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ABSTRACT 


This  document  describes  the  energy  distribution  within  an  electrically  small  antenna.  Chu 
(1948)  derived  expressions  for  the  “Q”  of  an  omnidirectional  antenna  whose  entire  structure 
is  enclosed  within  a  sphere.  Many  authors  have  extended  this  technique.  Because  the  sphere 
is  not  conformal  to  the  structure  of  most  linear  polarized  antennas,  much  energy  is  not 
accounted  for.  If  the  structure  is  surrounded  by  an  ellipsoid  of  revolution,  the  fit  is  much 
better.  Prolate  spheroidal  coordinates  are  expected  to  provide  a  much  better  description  for  a 
linearly  polarized  antenna.  We  derive  expressions  for  the  electromagnetic  field  and  the 
Poynting  vector  for  vertically  and  horizontally  polarized  antennas.  From  the  expressions,  we 
calculate  the  Q  of  an  electrically  small  antenna.  We  can  derive  an  analytic  expression  for  the 
value  of  Q  for  the  leading  term.  The  leading  term  is  a  product  of  the  inverse  of  the  cube  of 
the  dimensionless  wave  number  (product  of  the  wave  number  and  half  the  distance  between 
foci  of  the  ellipse)  and  a  factor  that  depends  only  upon  the  shape.  This  result  is  reformulated 
so  that  as  the  shape  changes,  the  volume  and  wave  number  remain  constant. 
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INTRODUCTION 


Since  1 947,  many  papers  have  been  published  on  the  radiation  Q  of  electrically  small 
antennas.  With  the  exception  of  the  original  work  by  Wheeler  (1947),  most  of  these  papers 
involved  antennas  contained  within  a  sphere.  In  two  of  these  papers.  Watt  (1967)  and 
Simpson  (2004)  used  circuit  concepts  to  investigate  antennas  with  cylindrical  symmetry. 
Collin  and  Rothschild  (1964)  included  an  infinite  cylinder.  This  document  derives  a  result  for 
the  radiation  Q  of  electrically  small  antennas  contained  within  a  prolate  spheroid  to  under¬ 
stand  the  effect  of  aspect  ratio  defined  as  the  height  to  width. 


THE  INTERPRETATION  AND  IMPORTANCE  OF  RADIATION  Q 

As  in  the  quality  factor  for  a  circuit,  the  radiation  Q  of  an  antenna  is  defined  as  2n  times 
the  ratio  of  the  maximum  energy  stored  to  the  total  energy  radiated  per  cycle.  Radiation  Q 
implies  that  only  radiated  power  is  considered  dissipated.  The  radiation  Q  of  an  antenna  is  an 
important  parameter.  The  available  impedance  bandwidth,  radiation  efficiency,  and  power¬ 
handling  capability  are  inversely  proportional  to  radiation  Q.  Thus,  knowledge  of  the  mini-  ■ 
mum  radiation  Q  for  a  given  antenna  is  important  for  design.  The  impedance  bandwidth 
limitation  implied  by  radiation  Q  can  be  overcome  to  some  extent.  The  efficiency  and  power¬ 
handling  limitations,  however,  cannot  be  overcome. 

Geyi,  Jarmuszewskik,  and  Qi  (2000)  showed  that  a  one-port  antenna  is  equivalent  to  a 
lossy  network.  The  Foster  reactance  theorem  applies  and  Q~f/BW  provided  Q»l,  where  f  is 
the  resonant  frequency  and  BW  is  the  3-dB  bandwidth  of  the  impedance.  This  bandwidth 
assumes  a  single  reactive  element  tuning  the  antenna.  The  bandwidth  determined  in  this  way 
is  really  the  bandwidth-efficiency  product.  Efficiency  can  be  decreased  to  increase  bandwidth 
(Simpson,  2004).  Multiple  tuning  can  significantly  increase  the  impedance  bandwidth. 
DeSantis’  (1977)  theoretical  results  showed  that  physically  realizable  lossless  circuits  limited 
by  Foster’s  Reactance  Theorem  increase  the  impedance  bandwidth  of  an  electrically  small 
monopole  by  more  than  three  orders  of  magnitude. 

The  bandwidth  limitations  for  transmission  can  be  overcome  in  other  ways.  One  way  is  to 
change  the  coupling  to  the  antenna  (Wheeler,  1947).  A  recent  version  of  this  approach  used 
at  U.S.  Navy  low-frequency  (LF)  stations  uses  a  transmission  line  and  network  to  transform 
the  series  resonance  of  the  antenna  to  parallel  resonance  fed  by  a  low-impedance,  solid-state 
amplifier.  The  amount  of  bandwidth  that  can  be  enhanced  by  this  method  is  a  function  of  the 
excess  KVA  available  in  the  transmitter.  Another  method  of  bandwidth  enhancement  at  the 
U.S.  Navy’s  very  low-frequency  (VLF)  station  at  Cutler,  Maine,  uses  variable  reactive 
energy  storage  elements  (Simpson,  2004). 

The  relationship  between  the  radiation  Q  and  radiation  efficiency,  r),  is  given  by 

n  =  Rr/(Rr+R,)=l/(l+Qr/Qt),  (1) 

where  Rr  is  the  radiation  resistance,  Ri  is  the  loss  resistance  in  the  tuning  element  (and  other 
losses),  Qr  is  the  radiation  Q,  and  Qt  is  the  Q  of  the  tuning  element  including  the  other  losses. 
Equation  (1)  asserts  that  a  lower  radiation  Q  implies  greater  radiation  efficiency.  A  higher  Q 
implies  higher  voltages  and  currents  as  well  as  higher  reactive  fields  near  the  antenna.  The 
higher  reactive  fields  result  in  increased  loss. 
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For  a  given  antenna,  the  voltage-limited  radiated  power  is  directly  related  to  the  radiation 
Q.  The  maximum  voltage  on  the  antenna  is  proportional  to  the  product  of  Qr  and  the  square 
root  of  the  power  radiated.  For  a  given  power,  this  voltage  depends  only  on  the  radiation  Q, 
not  the  loaded  Q.  Lower  Q  is  desired. 

PREVIOUS  RESULTS 

Wheeler  (1947)  was  the  first  to  investigate  the  concept  of  the  minimum  radiation  Q  that 
can  be  obtained  for  an  electrically  small  antenna.  He  introduced  the  term  “radiation  power 
factor”  as  the  inverse  of  the  radiation  Q.  Using  circuit  concepts  for  a  capacitor  and  inductor 
acting  as  an  antenna,  he  showed  that  the  radiation  power  factor  for  either  kind  of  antenna  is 
somewhat  greater  than 

Qr1  >  4jt2/3  •  Vol/X3,  (2) 

where  Vol  is  the  volume  of  the  cylinder  containing  the  antenna  and  X  is  the  wavelength.  The 
same  circuit  concepts  can  be  extended  to  show  that  the  voltage-limited  radiated  power  of  an 
electrically  short,  top-loaded  monopole  is  proportional  to  the  square  of  the  area  (Watt,  1967). 

Chu  (1948)  was  the  first  to  derive  analytic  results  using  Maxwell’s  equations  for 
omnidirectional,  linearly  polarized  radiation.  He  obtained  expressions  for  the  minimum 
radiation  Q  that  can  be  obtained  for  an  antenna  that  fit  within  a  sphere  of  a  given  radius.  This 
result  assumed  no  energy  remained  within  the  sphere.  The  fields  outside  the  sphere  were 
expanded  in  spherical  wave  functions.  These  modes  were  orthogonal.  The  stored  energy  and 
radiated  power  were  calculated  for  each  mode  using  the  complex  Poynting  Theorem  and  a 
ladder  circuit  derived  for  each  mode.  Chu  obtained  results  for  transverse-magnetic  (TM)  and 
transverse-electric  (TE)  linearly  polarized  and  circularly  polarized  modes.  Harrington  (1960) 
extended  this  theoretical  result  to  include  the  maximum  possible  gain.  Other  authors  have 
expanded  on  this  approach  and  derived  exact  results  (Fante,  1969;  McLean,  1996;  Grimes 
and  Grimes,  1995;  and  Fante  and  Mayhan,  1968)  and  have  used  the  same  type  of  theory  to 
examine  the  power-handling  capability  of  electrically  small  antennas  by  minimizing  the 
electric  field  on  a  sphere  around  the  antenna. 

The  literature  is  still  somewhat  controversial  on  the  Q  of  electrically  small  antennas, 
especially  for  multiport  antennas  (Grimes  and  Grimes,  2001).  For  electrically  small  single¬ 
port  antennas  contained  within  a  sphere,  the  results  agree.  When  the  sphere  radius  is  less  than 
one-tenth  of  a  wavelength,  all  results  reduce  to  the  minimum  possible  radiation  Q  given  by 

Qrmin'1  =  (ka)3  =  6JI2  Vol/X3,  (3) 

where  k  is  the  wavenumber  (2n/X),  a  is  the  sphere  radius,  and  Vol  is  the  volume  of  the  sphere. 

Note  that  equation  (3)  is  given  for  the  minimum  Q.  This  quantity  is  less  than  that  derived 
by  Wheeler  (1947)  in  equation  (2).  Since  Wheeler’s  derivation  was  based  on  circuit  theory 
considerations,  and  included  energy  stored  within  the  cylindrical  volume,  this  result  is 
reasonable. 

The  minimum  Q  for  equal  TM  and  TE  modes  phased  to  give  circular  polarization  is  half  of 
that  given  in  equation  (3).  Grimes  and  Grimes  (2002)  have  derived  results  for  antennas 
radiating  a  combination  of  TE  and  TM  modes.  They  showed  that  the  minimum  Q  depends  on 
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the  relative  phase  of  the  driven  elements.  They  gave  an  example  of  a  turnstile  antenna  that 
with  appropriate  phasing  has  a  minimum  Q  equal  to  one-third  that  given  in  equation  (3). 

Most  antennas  are  better  fit  by  cylindrical  enclosure  rather  than  a  sphere.  This  shape  gives 
two  variables,  height  and  diameter.  Equations  (1)  and  (3)  do  not  give  any  insight  into  the 
relative  effect  of  height  and  diameter  for  this  case.  Laport  (1952)  discussed  this  issue  with  the 
implication  that  VLF  antennas  should  be  built  as  high  as  possible  and  then  as  wide  as  neces¬ 
sary.  The  analysis  below  uses  prolate  spheroidal  coordinates,  which  somewhat  quantifies  this 
relationship. 

The  approach  below  is  similar  to  that  used  by  the  previous  authors.  The  fields  are  expand¬ 
ed  in  prolate  spheroidal  coordinates  and  the  energy  stored  and  the  power  radiated  are 
calculated  for  the  region  outside  an  ellipsoid  of  revolution  of  a  given  size.  We  show  that 
Maxwell’s  equations  are  separable  in  this  coordinate  system.  We  believe  this  is  a  new  result. 
Unfortunately,  the  modes  are  not  orthogonal  in  terms  of  power.  The  results  are  not  as  general 
as  the  previous  authors.  The  results  are  valid  for  the  minimum  radiation  Q  when  only  one 
mode  is  present.  We  derive  an  analytic  result  for  the  leading  term  in  the  expression  for  the 
minimum  radiation  Q  of  the  lowest  order  mode. 


GEOMETRICAL  RESULTS  FOR  PROLATE 
SPHEROIDAL  COORDINATES 

The  relation  between  Cartesian  and  prolate  spheroidal  coordinates  is  given  by 

x  =  0.5  d  cos((p)  V(u2  -  1)  V(l-r|2)  (4) 

y  =  0.5  d  sin((p)  V(  u2  -  1 )  V(  1  -r|2)  (5) 

z  =  0.5  d  u  r|  (6) 

This  document  includes  several  departures  in  nomenclature  from  Flammer  (1957).  The 
independent  radial  coordinate,  £  in  Flammer’ s  book,  is  represented  by  u  in  this  document. 

The  eigenvalue,  X  in  Flammer’ s  book,  for  the  angular  functions,  is  represented  by  k.  The 
variable  name,  X,  is  reserved  for  the  wavelength.  The  parameter,  7td/  X,  is  p  rather  than  c. 

The  variable,  c,  is  reserved  for  the  speed  of  light  in  vacuum.  In  one  equation  in  the  previous 
section,  the  variable,  q,  was  used  for  the  efficiency  of  the  antenna.  In  all  subsequent  sections, 
this  parameter  is  the  independent  variable  equal  to  the  cosine  of  the  polar  angle.  The  parame¬ 
ter,  d,  gives  the  distance  between  the  foci  of  the  ellipse.  The  equation  for  an  ellipse  with 
semi-minor  axis,  a,  and  semi-major  axis,  b,  is  given  by 

(x/a)2  +  (z/b)2  =  1.  (7) 

Substituting  the  expressions  for  x  and  z  from  equations  (4)  and  (6)  (with  cp,  the  azimuth 
coordinate  set  to  0)  into  equation  (7)  shows  a  connection  between  the  minimum  value  of  u 
(um>l)  and  the  ratio  of  b  to  a,  the  aspect  ratio.  The  surface  of  the  ellipsoid  of  revolution  is 
defined  by  u  =  um.  The  coordinate,  u,  varies  from  um  to  oo.  The  relation  between  um  and  the 
ratio  of  b  to  a  is  given  by 

um  =  (b/a)/V(b2/a2  -  1).  (8) 
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(9) 


The  distance  between  foci  of  the  ellipse,  d,  is  given  by 

d  =  2  b/  um  =  2  aV(b2/a2  -  1). 

The  expressions  can  be  simplified  slightly  by  defining  the  eccentricity  of  the  ellipse. 

An  antenna  with  length,  L,  and  radius,  r,  can  be  enclosed  either  by  a  sphere  or  by  an 
ellipsoid  of  revolution.  The  volume  of  the  antenna  is  given  by  kLt2.  The  minimum  volume  of 
the  sphere  that  surrounds  the  cylinder  is  given  by 

Vol  (sphere)  =  (4tt/3)((L/2)2  +  r2)3/2  ~  ttL3/6  +  7ir2L.  (10) 


in  which  the  approximation  of  a  long,  thin  cylinder  (L/2»r)  has  been  applied.  The  “wasted” 
volume,  i.e.,  the  volume  of  the  sphere  outside  the  cylinder  is  given  by  7tL3/6.  The  volume  of 
ellipsoid  of  revolution  that  encloses  the  cylinder  is  given  by 

Vol  (ellipsoid)  =  (4tc/3)  r2  b  /(I -(0.5  •  L/b)2).  (1 1) 


The  value  of  b  that  minimizes  the  volume  is  given  by  V3  L/2.  The  minimum  volume  is 
given  by 

V ol  (minimum  for  ellipsoid)  =  ir  r2  L  V3 .  (12) 


If  r  is  much  less  than  L/2,  this  minimum  volume  is  much  less  than  that  for  a  sphere.  The 
ellipsoid  of  revolution  fits  the  long,  thin  cylinder  much  more  closely  than  a  sphere.  In  the 
electromagnetic  case,  the  energy  within  the  ellipsoid  should  be  much  less  than  that  in  the 
sphere.  This  reduced  energy  implies  that  the  Q  derived  in  prolate  spheroidal  coordinates 
should  be  less  than  that  for  the  sphere.  Figure  1  shows  an  example  of  the  geometry.  The 
figure  of  the  ellipse  is  chosen  with  parameters  to  minimize  the  volume  of  the  ellipsoid  of 
revolution  that  encompasses  the  cylinder  representing  the  dipole.  The  values  for  the  parame¬ 
ters  of  the  cylinder  were  chosen  to  accentuate  the  idea  that  the  sphere  and  ellipsoid  of 
revolution  encompass  the  dipole. 
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The  parameter,  d,  is  the  distance  between  foci  of  the  ellipse.  Spherical  coordinates  are 
recovered  by  taking  the  simultaneous  limits,  d  — >  0  and  u— >x>,  while  the  product,  u  d  /2, 
approaches  r.  Similarly,  the  radius,  a,  is  approached  as  the  product,  umd/2,  in  which  um 
becomes  large  and  d  becomes  small. 


FIELDS  EXPRESSED  IN  PROLATE  SPHEROIDAL 
COORDINATES— VERTICAL  POLARIZATION 


Maxwell’s  equations  in  prolate  spheroidal  coordinates  are  most  conveniently  solved  by 
using  the  formalism  of  the  Hertz  potential.  In  this  method,  the  fields  are  the  second  deriva¬ 
tives  of  a  vector  function.  The  fields  are  the  first  derivative  of  the  usual  vector  and  scalar 
potentials.  The  addition  of  a  second  derivative  provides  a  freedom  that  allows  a  derivation  of 
Maxwell’s  equations  as  the  third  derivative  of  one  scalar  function.  The  formalism  of  Cartan 
frames  and  Debye  potentials  can  also  be  used  to  derive  the  results  in  a  coordinate  independ¬ 
ent  way  (Cohen  and  Kegeles,  1974).  This  formalism,  though  extremely  powerful,  provides  no 
additional  insight  for  analyzing  Maxwell’s  equations  in  the  Euclidian  spaces  appropriate  for 
ordinary  conditions. 

We  express  the  measured  fields,  E  and  B,  in  the  usual  way,  in  terms  of  the  vector  and 
scalar  potentials,  A  and  <E>,  by 

E  = -At  -  grad  (d>)  (13) 

B  =  curl  (A)  (14) 


A  bold  letter  denotes  a  vector  and  a  subscript  denotes  differentiation  with  respect  to  the 
variable.  The  vector  and  scalar  potentials  are  related  to  the  Hertz  potential,  II  by  the  relation: 

A  =  n,/c2  (15) 

0>  =  -div(n)  (16) 


By  expressing  the  vector  and  scalar  potentials  in  this  way,  the  Lorentz  gauge  is  satisfied 
identically,  and  a  wave  equation  results.  The  symbol,  c,  denotes  the  speed  of  light  in  vacuum. 

In  any  curvilinear  coordinate  system,  the  divergence,  gradient,  and  curl  are  related  directly 
to  the  metric,  ds2,  unit  vectors  u,  a,  (£,  and  the  metric  tensor,  hi,  which  has  only  diagonal 
coefficients  because  of  the  orthogonality  of  the  coordinate  system,  by  the  relation, 


dy2  +  dz2  =  hj2  du2  + 1122  dq2  +  h32d(p2 

(17) 

hi  =  0.5  d  V(u2-q2)/V(u2-l) 

(18a) 

h2=  0.5  d  V(u2-q2)/V(l-  q2) 

(18b) 

h3  =  0.5  dV(u2-l)  V(l-q2) 

(18c) 

The  independent  coordinate,  u,  is  defined  for  all  u  greater  than  um  >  1 .  The  coordinates,  q 
and  tp,  are  defined  for  all  values  between  -1  to  1  and  0  to  2n,  respectively.  We  will  call  u  the 
radial  coordinate  and  q  the  angular  one. 
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The  gradient  is  given  by 

grad(O)  =  u  d>u  /h,  +  n  <JVh2  +  (£  <t>,p  /h3  ( 1 9) 

The  divergence  is  given  by 

div  (A)  =  ((h2  h3  u-A)u  +  (h,  h3  tpA),,  +  (h,  h2  (fc-A),,  )/(h,  h2  h3)  (20) 

The  curl  is  given  by 

curl  (A)  =  u  ((h3  (g-A^  -  (h2  a-A),,,)/^  h3)  +  a  ((hi  U'A),  -  (h3  <£-A)u)/(h3  hi)  + 

< 2  ((h2  a-A)u  -  (hi  u-A)n  /(h,  h2)  (21) 

Specialization  to  a  vertically  polarized  system  occurs  if  (Flammer,  1957,  p.  70): 

A  =  (u  r|  /hi  +  a  u  /  h2)  0.5  d  et  /c2  (22a) 

O  =  -0.5ddiv((ur|/hi +au/h2)e)  (22b) 


The  choice  of  the  functional  form  for  A  and  d>  in  terms  of  u  and  q  comes  from  the  defini¬ 
tion  of  a  unit  vector  in  the  z  direction,  (z).  The  vector  potential  should  be  proportional  to  z 
so  that  the  electric  field  also  has  this  polarization.  We  also  specialize  to  omnidirectional 
antennas  by  requiring  that  no  quantity  depends  upon  the  azimuth  coordinate,  (p.  The  depend¬ 
ent  variable,  e,  is  a  scalar  quantity  and  an  arbitrary  function  of  t,  u,  and  q.  The  current  densi¬ 
ty,  J,  in  Maxwell’s  equations  must  also  be  in  the  vertical  direction  so  that  we  also  have 

J  =  (u  q  /hi  +  £  u  /  h2)  0.5  d  vt  /c2  (22c) 

This  equation  introduces  the  function  v  and  its  time  derivative.  Inputting  the  assumptions 
presented  in  equations  (22a),  (22b),  and  (22c)  into  the  definitions  of  the  gradient,  divergence, 
and  curl  leads  to  two  identical  equations  for  the  time  and  space  variation  of  e.  After  integra¬ 
tion  with  respect  to  t,  the  one  independent  equation  is 

(p0  V  —  e,t  /  c2)  d2  (u2-r)2)/4  =  -(u2- 1 )  euu  -  ( 1  -q2)  c,)T]  +  2qen  -  2  ueu  (23) 

Both  components  of  the  equation  for  e  result  in  the  same  equation.  If  a  time  dependence 
for  e  is  given  by  exp(-j  co  t)  and  the  current  density  is  zero,  the  solution  to  equation  (23)  is 
given  by  the  product  of  eigenfunctions  of  spheroidal  wave  equations  for  m  =  0: 

0  =  (( 1  -q2)S(n,0)ri)n  +  (K(n’0)  -  p2  q2)S(n’0)  (24a) 

0  =  ((u2- 1  )R(n’0)u)u  +  (p2  u2  -  K(n’0))R(n’0)  (24b) 

The  separation  constant,  K(n  0),  is  the  eigenvalue  of  equation  (24a)  adjusted  so  that  the 
solution  is  regular  for  all  q  between  -1  and  1 .  The  eigenvalue  is  a  function  of  the  frequency, 
co,  and  its  related  parameter,  p  =  co  d/(2c).  The  Appendix  describes  the  numerical  solution  of 
equations  (24a)  and  (24b).  Tables  for  the  eigenvalues  and  normalized  eigenfunctions  are 
described  in  Flammer  (1957)  and  Abramowitz  and  Stegun  (1965).  In  all  subsequent  work,  we 
assume  simple  exponential  time  dependence.  Figure  2  compares  the  numerical  results  for 
solving  the  eigenvalue  equation  with  those  in  Flammer  (1957).  The  comparison  always 
agrees  with  every  digit  listed  in  Flammer’s  book.  The  Appendix  presents  the  methods  used  in 
computing  the  eigenvalue,  eigenfunctions,  or  the  angular  and  radial  equations. 
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Figure  2.  Eigenvalue  computation. 


In  terms  of  e,  the  electric  field,  E,  and  the  magnetic  field,  B,  are  given  by: 

E  =  T1  u  {P2  e  +  ((e(u2-l)V(u2-ii2))u  +  (u  (e(l-if))n/(u2-if))u  /r, } (4/d2) V(u2- 1 )/ V(u2-rj2)+ 
u  H  {P2  e  +  0i(e(u2- 1  ))u/(u2-T|2))r|  /u+  ((e(  1  -r]2))T/(u2-ri2))n }  (4/d2)V(  1  -rj2)/V(u2-rj2)  (25a) 

B  =  -4  j  p  (2  (ueu-iie,)  •  {cd2(u2-if)}-V(u2-l)  V(l-r|2)  (25b) 

These  expressions  are  used  in  the  subsequent  section  to  calculate  the  Poynting  vector.  We 
will  eventually  expand  the  dependent  variable,  e,  in  a  series  of  the  complete  set  of  eigenfunc¬ 
tions.  The  radial  function,  R,  will  be  specialized  to  deal  only  with  outgoing  waves. 


FIELDS  EXPRESSED  IN  PROLATE  SPHEROIDAL 
COORDINATES— HORIZONTAL  POLARIZATION 

Calculation  of  the  fields  and  the  Poynting  vector  are  much  simpler  for  the  case  of 
.  horizontal  polarization  than  for  vertical  polarization.  This  calculation  is  important  if  circular 
polarization  is  required.  The  Hertz  potential  has  only  one  component, 

A  =  <2  et  /  c2  (26a) 

<S  =  0  (26b) 

J  =  <2  vt  (26c) 

Applying  the  formulas  leads  to  one  non-trivial  equation: 

-Po  V  d2  (u2  -  rf)/4  =  p2e  (u2  -  if)  +  (u2  -  1  )euU  +  2ueu  -  e/(u2  -  1 )  +  ( 1  -  q2)en^  -  2r|en  —  e/(  1  -  if) 

(27) 
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The  solution  to  equation  (27)  leads  to  the  set  of  eigenfunctions  for  the  spheroidal  wave 
equations  for  m  =  1 : 

0  =  (( 1  — iq^)S(n’ 1)T|)ij)  +  (K(n'u  -  p2  q2  -  1/(1  -  r|2))  S(n>1)  (28a) 

0  =  ((u2-l)R(n,1)u)u)  +  (P2  u2  -  K(n,1)  -  l/(u2  -  l))R(n,1)  (28b) 

As  in  the  vertical  polarization  case,  only  outgoing  waves  are  chosen.  In  the  language  of 
prolate  spheroidal  coordinates,  the  radial  function  will  be  asymptotic  to 

R(n,l)  ~  (1/(P  u))  exp(j  (P  u-(n  +  1)ji/2)  as  u— >oo  (29) 

This  expression  is  called  the  radial  function  of  the  third  kind.  The  complex  conjugate  of  R 
is  called  the  radial  function  of  the  fourth  kind. 

The  expressions  for  the  electric  and  magnetic  fields  are  given  by 

E  =  (2p/d)2  (g  e  (30a) 

B  =  -4  j  (p/(cd2))  {u  ((l-q2)en  -  qe/(l-q2))V(l-q2)  V(u2-q2)  -  a  ((u2-l)eu  +  e/(u2-l))V(u2-l)/ 

V(u2-q2)}  (30b) 


EXPRESSIONS  FOR  POYNTING  VECTOR- 
VERTICAL  AND  HORIZONTAL  POLARIZATIONS 

Two  surfaces  are  important  for  this  problem.  The  first  is  the  one  that  surrounds  the 
antenna,  which  is  defined  by  u  =  um  >  1 .  The  second  is  located  at  distances  that  are  large 
when  compared  to  a  wavelength  (u»l).  Only  that  component  of  the  Poynting  vector  in  the  u 
direction  will  be  important  for  describing  radiation.  The  other  component  is  tangential  to 
these  surfaces. 

The  expression  for  the  Poynting  vector,  S,  the  power  per  unit  area  through  a  surface,  is 
given  by 

S  =  E  x  B*/(2p0)  (31) 

The  power  through  the  surface  defined  by  u  =  constant  is  given  by 

i 

P  (u)  =  2ti  j  h2  h3  u  (E  x  B*)/(2p0)  dr)  (32) 

-l 

Inserting  the  expressions  for  h2,  h3,  E,  and  B  results  in 

P  (u)  =  (2  7i  j  co  u  (u2- 1  )/p0dc2)j  dq  (1-q2)/  (u2-q2)2  (p2u2e  +  (u2-l)(u2+  rf)eu/(u(u2-  q2))+ 
r|eu11  (u2-l)/u+e^(l-q2)-2  qen-p2  q2  e+2  q  e^  ( 1  -  q2)/(u2-  q2)}  (ueu*-  qen*)  (33) 

In  equation  (33),  the  assumption  of  a  time  dependence  of  exp(-zcot)  has  been  made.  If  we 
expand  the  dependent  variable,  e,  in  a  series  of  the  complete  set,  S(n,0)  as 

OO 

e  =  £  en  S(n’0)(f3,q)  R(n,0)(P,u)  (34) 

n  =  0 
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equation  (33)  then  becomes 


P  (u)  =  (2  n  j  to  u  (u2-l)/p0dc2)  ££en  e,*|  dq  (l-q2)/(u2-q2)2  {p2u2S(n’0)R(n’0)-  KnS(n’0)  R(n'0)+ 

11=0, 1=0 

(u2-  1  )(u2  +  q2)/(u(u2  -  q2)S(n,0)R(n’0)u  +  q(u2  -  l)S(n’0)„R(n'0)u/u  +  2q(l  -  q2)R(n'0)S(n’°V 
(u2  -  q2)}  {uS(l’0)R(l,0)*u  -  qR(1’°>*S(W),}  (35) 

In  the  limit  as  u— >co,  the  radial  function  take  their  asymptotic  forms  (equation  (29))  and  the 
total  power  for  vertically  polarized  radiation  becomes 

P  (u  >oo)  =  (47tp2/d2cp0)  £ien  e,*J  dq  (1-q2)  S(n>0)S(1’0)  (36) 

n=0,l=0 


Only  if  only  one  mode  is  excited  is  the  total  power  positive  definite.  Unfortunately,  the 
orthogonality  of  the  S  in  the  weight  function  of  equation  (36)  cannot  be  assumed.  The 
presence  of  terms  with  u  in  the  denominator  makes  the  orthogonality  doubtful.The  situation 
for  horizontally  polarized  radiation  is  much  simpler.  Inserting  the  expressions  for  the  electro¬ 
magnetic  field  into  the  expressions  for  the  power  leads  to 

P  (u)  =  (27tjco3d/4c4p0)  Jdq  e  ((u2-l)  eu*  +  u  e*)  (37) 

Expanding  the  dependent  variable,  e,  in  a  series  of  S(n,1)  leads  to 

P  (u)  =  (47tjp3/d2cp0)X  |  e„  | 2  R(n,1){uR(n,1)*  +  (u2-l)R<nJ)*u)  Jdq  S(n’1)2  (38) 

n=l 

The  double  sum  reduces  to  a  single  one  because  of  the  orthogonality  of  the  S  over  the 
interval.  In  the  limit  of  large  u,  this  expression  for  horizontal  polarized  radiation  reduces  to 

Power  (u — »oo)  =  (4np2/d2cq0)  £  |  en  |  2Jdq  S(n>1)2  (39) 

n=l 

The  power  at  infinity  is  a  positive  definite  function.  The  expression  in  equation  (39)  goes 
directly  over  to  the  analysis  given  in  Chu  (1948)  and  Grimes  and  Grimes  (1995)  for  deter¬ 
mining  the  expression  for  the  Q  of  the  antenna. 


APPROXIMATION  OF  SMALL  ELECTRICAL  SIZE 

If  small  electrical  size  is  approximated,  some  formulas  can  be  simplified.  The  parameter 
that  relates  to  electrical  size  is  P  =  ©  d/(2c).  Small  p  is  equivalent  to  small  electrical  size. 
Most  of  the  formulas  are  functions  of  the  square  of  p. 

The  equation  for  the  angular  eigenfunction  in  prolate  spheroidal  coordinates  is 

0  =  ((l-q2)S(n’m\)„)  +  (K(n’m)  -  P2  q2  -  m2/(l-q2))  S(n’m)  (40) 

This  equation  is  almost  identical  to  that  for  associated  Legrendre  polynomials: 

0  =  ((l-q2)P„mn)„)  +  (n(n+l)  -  m2/(l -q2))  Pnm  (41) 

Only  the  term  p2  q2  is  different.  This  term  should  be  considered  a  small  correction  for  an 
electrically  small  antenna.  Because  of  the  standard  notation  for  Legendre  polynomials,  we 


9 


have  put  the  subscript,  1,  on  P|m.  This  should  not  cause  confusion  with  the  derivative  with 
respect  to  a  variable. 

Abramowitz  and  Stegun  (1965,  p.  334)  give  the  recursion  relation  for  the  Legendre 


polynomials: 

(n-m+1)  Pn+1m  =  (2n+l)  p  Pnm  -  (n+m)  Pn_,m  (42) 

Using  this  recursion  relation,  the  orthogonality  of  the  Legendre  polynomials,  and  the 
integral  of  the  square  (determined  to  the  normalization), 

1-1  dp  Pnm  Pnm  =  2  (n+m)!/((2n+l)  (n-m)!),  (43) 

we  determined  the  following  integrals: 
i 

Li  dp  p  Pnm  P„+im  =  2  (n+l+m)  (n+m)!  /((2n+l)  (2n+3)  (n-m)!)  (44a) 

J.i  dp  p  Pnm  Pn.im  =  2  (n-m)  (n+m)!/((2n-l)  (2n+l)  (n-m)!)  (44b) 

I.,'  dp  p2  Pnm  P„m  =  2(2n2+2n- 1  -2m2)(n+m)!/((2n+ 1 )(  2n-l)(  2n+3)(n-m)!)  (44c) 

Li'dp  p2Pnm  Pn+2m  =  2(n+2+m)(n+ 1  +m)(n+m) !/((2n+ 1  )(2n+3)(2n+5)(n-m) ! )  (44d) 

I 

J.i  dp  p2Pnm  Pn-2m  =  2(n-m)(n-l-m)(n+m)!/((2n-l)(2n-3)(2n+l)(n-m)!)  (44e) 


We  seek  to  determine  the  eigenfunctions  and  eigenvalues  of  equation  (40)  through  first 
order  in  p2.  Most  relations  are  well  determined  by  zeroth  order.  We  seek  to  understand  the 
effect  of  finite  but  small  size  upon  these  relations.  Accordingly,  we  express  the  eigenvalue, 


K(n,m),  as  a  two-term  expansion, 

K<nm)  =  n(n+l)  +  y(n,m)p2,  (45) 

and  the  eigenfiinction,  S(nm),  as  a  three-term  expansion, 

S(n’m)  =  do(n’m)  Pnm  +  P2  (d,(n’m)  Pn+2m  +  d2(n’m)  Pn-2m)  (46) 

Only  two  of  the  terms  contribute  in  equation  (46),  if  n  is  0  or  1 .  Substituting  equations  (45) 
and  (46)  into  equation  (40)  and  using  the  various  results  above  leads  to 

K(n'm)  =  n(n+ 1 )  +  ( 2n2+2n- 1  -2m2)  p2  /((2n- 1  )(2n+3))  (47) 

and 

d,(n’m)  =  -(n+2-m)(n+l-m)  do(n>m)  /(2(2n+3)2(2n+l))  (48a) 

d2(n>m)  =  (n+m)(n-l+m)  d0(n’m)  /(2(2n-l)2(2n+l))  (48b) 


The  overall  factor,  do(n’m),  is  determined  by  the  scale  factor  of  the  Legendre  polynomials 
(Flammer,  1957,  p.  21): 

S(n  m)  (P,p  =  0)  =  Pnm  (p  =  0)  =  (-l)((n  m)/2)(n  +  m)!/(2n((n  -  m)/2)!((n  +  m)/2)!)  (49a) 
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if  n-m  is  even,  or 

S(n,m)(P,r|  =  0) ,,  =  Pnm(ri  =  0)  „  =  (-l)((n'm‘l,/2)(n  +  m  +  l)!/(2n((n -  m  -  l)/2)!((n  +  m  +  l)/2)!), 

(49b) 

where  the  derivative  is  evaluated  at  r|  =  0  for  n  -  m  odd.  The  final  result  for  do<n,m)  is  compli¬ 
cated  and  will  not  be  written  here.  Because  the  eigenfunctions  are  expressed  as  three 
Legendre  polynomials,  they  are  easily  evaluated.  The  eigenvalue  is  a  simple  function  of  the 
parameters,  n,  m,  and  p. 

The  radial  function  is  somewhat  more  difficult  to  analyze.  In  general,  the  radial  function  is 
given  as  an  integral  over  q  of  a  weighted  product  of  spherical  Hankel  functions  and  the 
eigenfunctions  for  the  angular  equation.  The  general  expression  is  given  by 

R(n,ra)(P,u)  =  jn  m  (2m  +  l)(u2-  l)m/2/(2m+1rn!dmn0(p)) 

Li  'dr|h(1)m[P(u2  +  r|2  -  1),/2](1  -  t12)m/2S<n,m)(P,  r|)(u2  +  q2  -  l)"1”72  (50a) 

for  n  -  m  even,  or 

R(n,m)(P,u)  =jn'm‘'  (2m  +  3)u(u2  -  l)m  2/(2m  +  ,m!dmn,(p)) 

JYdq  r)  h(1)m+i[P(u2  +  ri2- 1)1/2](1  -  r12)m/2S(n’IT,)(p,  q)(u2+q2-  l)'(m+1)/2  (50b) 


for  n  -  m  odd. 

The  function  h(1)m  is  the  spherical  Hankel  function  of  the  first  kind.  In  the  limit  of  small 
argument,  the  Hankel  function  is  dominated  by  the  imaginary  part,  a  spherical  Neuman 
function.  For  small  values  of  the  argument,  the  Hankel  function  is  given  by  Abramowitz  and 
Steagun  (1965,  p.437): 

h(1)m  (z  )  =  -j-l-3-5-...(2m-l)  z(m+1)  [  l-(z2/2)/(l!(l-2m))+(z2/2)2/(2!(l-2m)(3-2m))...]  (51) 

In  the  limit  of  small  electrical  size,  the  parameters  dmn0  and  dmni  become  approximately 
equal  to  1  and  the  function  S(m,n)  becomes  an  associated  Legendre  polynomial.  The  integral  in 
equation  (50a)  is  easily  done  for  m  =  0  and  n  =  0.  The  result  is 

R(0  0)(P,  u)  =  -j/(2p)  ln((u+l)/(u-l))  (52) 

The  integral  can  also  be  evaluated  for  m  and  n  less  than  2  in  terms  of  elementary  functions. 
The  results  can  probably  be  extended  to  almost  any  order  by  using  differentiation  under  the 
integral  sign  and  the  following  identity: 

i  dq  (s+wq2)'172  =  w'172  ln(q  w172  +  (s+wq2)172)  (53) 

Differentiation  with  respect  to  s  in  equation  (53)  leads  to  higher  powers  of  (s+wq2)  in  the 
denominator.  Differentiation  with  respect  to  w  leads  to  higher  even  powers  of  q  in  the 
numerator  and  higher  powers  of  (s+wq2)  in  the  numerator. 

Comparisons  with  numerical  evaluation  of  the  integral,  equation  (50a),  and  solution  of  the 
original  differential  equation  (24b)  show  that  this  expression  and  that  for  the  derivative  with 
respect  to  u  are  accurate  to  three  digits  if  the  value  of  u  is  less  than  0.3/p.  Accuracy  for  larger 
u  is  extended  if  the  form  chosen  is  given  by 

R(0’0)(P,  u)  =  -0/2  P)  expO  P  u)  ln((u  +  l)/(u  -  1))  (54) 
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Figures  3,  4,  5,  6,  and  7  compare  the  numerical  evaluation  of  the  differential  equation  and 
the  analytic  function  in  equation  (53). 


Figure  3.  Numerical  and  analytic  approximation  comparison  for  R(3). 


Figure  4.  Error  in  computing  real  and  imaginary  parts  from  analytic  formula  (u  =  1 .005). 
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Figure  5.  Error  in  computing  real  and  imaginary  parts  from  analytic  formula  (u  =  1 .01). 


Figure  6.  Error  in  computing  real  and  imaginary  parts  from  analytic  formula  (u  =  1 .05). 
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Figure  7.  Error  in  computing  real  and  imaginary  parts  from  analytic  formula  (u  =  2). 


In  these  figures,  we  have  evaluated  the  relative  accuracy  while  keeping  u  constant  and 
varying  p  and  varying  u  while  keeping  p  constant.  This  form  for  the  radial  function  will  be 
used  in  subsequent  sections.  The  accuracy  is  sufficient  for  our  purposes. 


ANTENNA  Q 

The  Q  of  an  electronic  device  is  equal  to  the  ratio  of  bandwidth  of  operation  to  the  center 
frequency.  Antenna  bandwidth  has  multiple  definitions.  They  include  full  width  at  half 
maximum  of  the  gain  versus  frequency,  frequency  span  over  which  the  return  loss  is  10  dB 
lower  than  its  maximum,  or  the  span  over  which  the  voltage  standing  wave  ratio  (VSWR)  is 
lower  than  a  certain  number.  Each  definition  has  its  use.  We  will  use  the  relation  between 
power,  P(u),  radiated  through  a  surface,  and  the  energy,  W,  stored  within  the  volume  bound¬ 
ed  by  the  surface. 

The  power  has  already  been  evaluated  for  an  arbitrary  distribution  of  modes.  The  expres¬ 
sion  for  Q  is  given  by 

Q  =  2  (£>  Ws,ored/P(u— >Co)  (55) 

This  equation  is  a  common  expression  for  Q  when  dealing  with  antennas.  As  Chu  (1948) 
shows,  this  expression  can  be  related  to  the  parameters  of  a  circuit.  In  general,  the  expression 
for  the  energy  within  the  field  is  determined  from  Maxwell’s  equation  by 

Wtotai  =  JdV  (e0  E-EV4+  B-B7(4po))  (56) 

The  volume  element  in  prolate  spheroidal  coordinates  is  given  by 

JdV  =  Jhih2h3  dudr|d(p  =  2n{dJ2f\ (u2-q2)  du  dr|  (57) 
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The  integral  above  is  evaluated  over  all  space  external  to  the  volume.  The  idea  is  that  no 
energy  is  stored  close  to  the  antenna.  One  reason  for  using  prolate  spheroidal  coordinates  is 
that  there  is  less  excluded  volume. 

The  variable,  u,  is  between  um  and  oo,  and  q  is  between  -1  and  1 .  The  stored  energy  is  the 
difference  between  the  total  energy  given  in  equation  (56)  and  that  radiated  away.  The  energy 
per  unit  volume  radiated  away  is  the  Poynting  vector  divided  by  the  speed  of  light.  At 
infinity,  subtracting  the  value  of  the  Poynting  vector  cancels  the  asymptotic  value  of  the  total 
energy  so  that  the  total  stored  energy  is  a  finite  value. 

The  stored  energy  can  be  written  as 

o°  1 

W stored  =  2tc  (d/2)3  Jum  J  ,(u2  - q2)  du  dq(e0  EE*  /4  +  B*B7(4po)-u*E  x  B*/(2p0  c))  (58) 

The  expression  in  equation  (58)  is  a  variation  of  the  expression  in  Collin  and  Rothschild 
(1964)  or  McLean  (1996)  for  an  antenna  enclosed  by  a  sphere  of  radius,  a.  Using  the  expres¬ 


sions  for  E  and  B  already  given  in  previous  sections,  we  find 

W stored  =  n  e0  d  1  |um  duj  4  drj(rj2(u2- 1 )  ||Qi|2/(u2  -  ti2)4-Hu2(l  -  r|2)|Q2|2/(u2  -  q2)4  + 

p2(u2  -  1)  (1  -  q2)  IIQ3IV-TI2)2  -2j Q2Q3 *u(  1  -q2)  Vu2- 1  / (u2-q2)772  ,  (59) 

where 

Qi  =  R  S  (K-p2q2)+(l-  q2)Sr/  q(u  Ru-(u2+  q2)R/(u2-  q2))-2u(u2-l)SRu/(u2-  q2)  (60a) 

Q2  =  (P2u2-k)R  S  +(u2-l)Ru/u(qSn+(u2+q2)S/(u2-q2))+2qR(l-q2)Sn/(u2-q2)  (60b) 

and 

Q3  =  S  u  Ru  -  R  q  Sn  (60c) 

In  terms  of  the  dependent  variables  given  above,  the  power  radiating  through  the  surface  at 
u  =  um  is  given  by 

P  =  4n  j  p  u  (u2- 1  )/(d2cp0)J_]  'dq  (1  -q2)Q2Q37(u2-q2)3  (6 1 ) 


In  general,  the  integrals  must  be  done  numerically,  even  when  the  radial  and  angular 
functions  are  given  in  terms  of  simple  analytical  functions. 

Fortunately,  the  integrals  can  be  done  in  equation  (59)  for  the  leading  term  in  p.  This  term 
will  lead  to  an  inverse  cube  relationship  between  Q  and  the  dimensionless  wave  number.  This 
relationship  is  important  in  its  own  right.  The  result  also  gives  confidence  in  the  numerical 
results  presented. 

From  equation  (54)  we  use  the  expression  for  the  first  derivative  of  the  radial  function: 

Ru  (u)  =  (,5-ln((u+l)/(u-l))  +  j/(p  (u2-l)))  expO  P  u)  (62) 

The  small  p  limit  of  the  eigenvalue  and  the  angular  eigenfunction  are  given  by: 

X  =  p2  /3  (63a) 

S  =  1  -  p2q2/6  (63b) 
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The  expressions  above  are  convenient.  Any  derivative  of  S  with  respect  to  the  angular 
coordinate  leads  to  an  extra  power  of  p2.  The  eigenvalue  is  also  proportional  to  this  extra 
factor.  The  eigenfunction  can  usually  be  set  equal  to  1  with  little  error.  The  leading  term 
occurs  from  the  factor  of  P  in  the  denominator  of  the  expression  for  Ru.  The  similar  term  in 
the  denominator  of  R  in  equation  (54)  is  always  multiplied  by  p2. 

The  expression  for  the  stored  energy  from  equation  (59)  is  given  by 

Wstored  =  (7te0/dp2)Ium00duJ.i  'dr|  [  4  r|2u2(u2-  1)  +  (1  -  r|2)(u2  +  r|2)2](u2  -  r|2)  4  (64) 

The  integrals  are  done  in  a  straightforward  manner.  The  book  of  integrals  given  in  Dwight 
(1961,  pp.  35-36)  is  one  example.  The  final  result  is 

Wstorcd  =  (jie0/ dp2)  [um/(um2-  1)  -  0.5/um+  0.25(um2-  l)ln((um+  l)/(um-  1))]  (65) 

In  the  same  limit  the  total  power  is  given  by 

P  total  =  167tp2/(3d2cp0)  (66) 

Equation  (65)  represents  one  of  the  more  important  results  of  this  document.  Because  Q  is 
the  product  of  the  radial  frequency,  co  =  2cp/d  and  Wst0red/Ptotai,  the  final  result  for  Q  is 

Q  =  3/(4p3)  [um/(um2  -  1)  -  0.5/um  +  0.25(um  2-l)ln((um  +  l)/(um -  1 ))]  (67)' 

Figure  8  shows  the  result  for  Q  given  for  the  lowest  mode  as  a  function  of  P  for  six  values 
of  the  minimum  um.  For  values  of  p  less  than  0.01 ,  the  expression  for  Q  given  by  equation 
(67)  is  confirmed  by  doing  the  integral  in  equation  (64)  numerically  to  four  digits. 


Figure  8.  Logarithm  of  Q  versus  Beta. 


Table  1  lists  the  logarithm  to  base  10  of  the  value  of  Q  versus  the  dimensionless  wave 
number  for  six  values  of  the  shape  factor,  um.  If  this  shape  factor  gets  close  to  1,  the  shape 
becomes  more  like  a  line.  As  this  factor  increases,  the  shape  becomes  more  spherical. 
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Table  1 .  Logarithm  to  the  base  10  of  Q  versus  p  for  six  values  of  um. 


p\um 

1.005 

1.010 

1.050 

1.200 

1.500 

2.000 

0.010 

7.8761 

7.5718 

6.8610 

6.2029 

5.6834 

5.1989 

0.025 

6.6822 

6.3779 

5.6671 

5.0093 

4.4902 

4.0065 

0.040 

6.0697 

5.7654 

5.0547 

4.3972 

3.8788 

3.3966 

0.055 

5.6546 

5.3504 

4.6398 

3.9827 

3.4654 

2.9854 

0.070 

5.3401 

5.0359 

4.3255 

3.6691 

3.1532 

2.6759 

0.085 

5.0869 

4.7827 

4.0725 

3.4169 

2.9027 

2.4288 

0.100 

4.8747 

4.5706 

3.8606 

3.2060 

2.6938 

2.2240 

0.115 

4.6922 

4.3880 

3.6784 

3.0249 

2.5150 

2.0497 

0.130 

4.5320 

4.2278 

3.5186 

2.8663 

2.3591 

1.8988 

0.145 

4.3891 

4.0850 

3.3761 

2.7253 

2.2210 

1.7661 

0.160 

4.2602 

3.9562 

3.2477 

2.5985 

2.0973 

1.6483 

0.175 

4.1428 

3.8388 

3.1308 

2.4832 

1.9855 

1.5427 

0.190 

4.0349 

3.7310 

3.0235 

2.3778 

1.8837 

1 .4474 

0.205 

3.9350 

3.6312 

2.9243 

2.2806 

1.7905 

1 .3609 

0.220 

3.8422 

3.5384 

2.8321 

2.1905 

1 .7045 

1.2819 

0.235 

3.7553 

3.4516 

2.7459 

2.1066 

1.6250 

1 .2095 

0.250 

3.6737 

3.3700 

2.6651 

2.0282 

1.5511 

1.1428 

0.265 

3.5967 

3.2931 

2.5889 

1 .9546 

1.4822 

1.0812 

0.280 

3.5238 

3.2204 

2.5169 

1 .8853 

1.4178 

1 .0241 

0.295 

3.4546 

3.1513 

2.4487 

1.8198 

1 .3573 

0.9710 
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DISCUSSION 


The  value  of  um  is  related  to  the  aspect  ratio  of  the  ellipsoid  of  revolution.  This  relationship 
is  displayed  in  equation  (8).  This  number  is  independent  of  the  volume  and  relates  only  to  the 
geometry  of  the  basic  ellipse.  To  relate  the  Q  to  the  total  volume,  V,  and  to  the  wave  number, 
k,  of  the  radiation,  is  straightforward.  The  total  volume  of  an  ellipsoid  of  revolution  with 
major  axis  length,  2b,  and  minor  axis  length,  2a,  is  given  by 

Vol  =  47ia2b/3  =  47rd3(um2-l)um/24  (68a) 

Becaue  a  factor  of  p3  =  (k  d/2)3  is  in  the  denominator  of  the  expression  for  Q,  the  volume 
can  be  extracted  from  the  final  result.  This  extraction  will  allow  the  shape-dependent  effects 
to  be  separated  from  the  volume.  The  result  is  as  follows: 

Q  =  (7i/2k3Vol)  um(um2-l)  [um/(um2-l)  -0 . 5/um+0 .25 (um2- 1  )ln((um+ 1 )/ (um- 1 ))]  (68b) 

As  shown  in  Figure  8  and  equation  (67),  for  a  given  value  of  P,  the  Q  is  a  monotonically 
decreasing  function  of  the  shape  factor,  um.  This  decrease  is  largely  due  to  the  change  in  the 
volume.  As  um  approaches  1 ,  the  value  of  Q  increases.  This  increase  is  reasonable.  A  long, 
thin  dipole  has  a  higher  Q  than  a  broad  one.  The  broad  dipole  has  many  paths  in  which  to 
complete  the  circuit.  The  resonance  width  is  broadened  with  increasing  lateral  extent.  The 
value  of  Q  is  a  decreasing  function  of  p. 

If  the  shape  and  the  wavenumber  remain  constant  and  the  value  of  the  volume  is  decreased, 
the  value  of  Q  increases.  For  constant  volume  as  the  value  of  um  increases,  the  value  of  Q 
also  increases.  The  volume  is  becoming  more  spherical.  Storage  of  energy  in  the  volume 
increases.  For  constant  volume  as  um  approaches  1  (a  very  thin  dipole),  the  value  of  Q 
approaches  zero. 

The  inverse  cube  dependence  of  Q  on  the  wavenumber  is  typical  for  an  electrically  small 
antenna.  Many  authors,  starting  with  Chu  (1948),  have  calculated  this  dependence  as  the  first 
term  in  a  power  series.  Grimes  and  Grimes  (1995)  have  provided  a  more  accurate  calculation 
that  reduces  to  such  a  power  series  in  the  appropriate  limit.  The  coefficient  multiplying  the 
k'3  is  somewhat  smaller.  For  a  small  dipole  surrounded  by  a  sphere  of  radius,  a,  the  value  of 
Q  (McLean,  1996,  among  many  others)  is 

Q sphere  =  (k  a)’3  +  (ka)'1  (69) 

For  um  =  2,  the  first  term  for  an  ellipsoid  of  revolution  is  given  by 

Qempsoid  =  0.158-(2/kd)3  (70) 

The  coefficient  of  the  k‘3  is  much  smaller  for  the  ellipsoid  with  a  particular  aspect  ratio 
than  for  the  dipole  surrounded  by  a  sphere.  The  sphere  surrounds  a  much  larger  volume.  This 
“wasted”  volume  contributes  significantly  to  the  Q  in  this  theoretical  limit. 

Figure  9  compares  Chu’s  results  for  the  sphere  (which  has  an  aspect  ratio  of  1)  with  those 
of  the  prolate  spheroid.  The  leading  terms  of  the  sphere  and  the  ellipsoid  of  revolution  are 
proportional  to  the  cube  of  the  wavelength  and  the  inverse  of  the  volume.  Multiplying  the 
value  of  Q  by  the  product  of  the  volume  divided  by  the  cube  of  the  wavelength  allows  a 
presentation  that  depends  only  upon  shape.  The  value  for  the  ellipsoid  of  revolution  becomes 
that  for  the  sphere  as  the  aspect  ratio  gets  close  to  1 .  If  the  volume  is  held  constant,  the  value 
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of  Q  for  the  ellipsoid  is  always  less  than  that  for  the  sphere.  As  the  aspect  ratio  gets  larger 
and  the  antenna  gets  taller  to  keep  the  volume  the  same,  the  value  of  Q  approaches  %  of  that 
for  the  sphere. 


Aspect  Ratio  of  Ellipsoid 


Figure  9.  Ratio  of  Q  for  ellipsoid  to  that  of  sphere. 

If  the  height  of  the  antenna  remains  constant,  the  value  of  Q  for  the  prolate  spheroid  is 
always  larger  than  that  for  the  sphere.  This  value  is  in  accord  with  experience.  Making  the 
antenna  thinner  makes  the  system  more  resonant,  and  the  Q  becomes  larger.  We  substitute 
the  expression  in  equation  (9)  relating  the  height  of  the  ellipsoid,  b,  to  the  distance  between 
foci,  d,  in  the  expression  for  (3,  to  obtain 

Qe=  0.75  QJ2-Tvbf  um3  [um/(um2-l)  -0.5/um+0.25(um  2-l)ln((um+l)/(um-l))]  (71) 

We  divide  Qe  by  Chu’s  expression  for  the  sphere  and  set  b  =  r  to  obtain 

Qe/Qsphere  =  0.75  um3  [um/(um2-l)  -0.5/um+0.25(um  2-l)ln((um+l)/(um-l))]  (72) 

Figure  10  shows  the  formula  above  as  a  function  of  the  aspect  ratio.  As  um  becomes  large, 
the  expression  in  brackets  approaches  4/(3um3).  As  um  approaches  1,  the  aspect  ratio  becomes 
large  and  the  first  term  dominates. 
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Figure  10.  Q  for  ellipsoid  compared  to  sphere. 


In  deploying  an  electrically  small  antenna,  space  and  height  will  probably  be  constrained. 
For  the  reasons  given  in  the  introduction,  the  value  of  Q  should  usually  be  minimized.  These 
results  indicate  that  keeping  the  volume  constant  and  allowing  the  height  to  increase  leads  to 
a  25%  reduction  in  the  value  of  Q.  Keeping  the  height  constant  while  reducing  the  volume 
leads  to  a  large  increase  in  Q. 


CONCLUSIONS 

Treating  the  electromagnetic  fields  in  prolate  spheroidal  coordinates  permits  a  separation 
of  the  volume  and  shape-dependent  effects  upon  the  value  of  Q.  The  shape  can  change  from 
spherical  to  cylindrical  in  appropriate  limits.  The  volume  can  remain  fixed  during  this  trans¬ 
formation. 
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APPENDIX:  NUMERICAL  METHODS  FOR  EVALUATING 
PROLATE  SPHEROIDAL  EIGENVALUES  AND  EIGENFUNCTIONS 


We  seek  to  develop  a  method  for  evaluating  the  eigenfunctions  of  the  angular  and  radial 
coordinate.  For  the  prolate  spheroidal  case,  the  equations  are  identical.  The  range  of  the 
independent  variable  is  the  primary  difference  between  the  two. 

The  book  on  numerical  methods  by  Press  et  al.  (1989,  p.  602)  suggested  a  relaxation  meth¬ 
od  for  solving  for  the  eigenvalue,  including  a  program  listing  in  FORTRAN.  We  chose  a 
somewhat  different  method,  one  that  would  provide  the  eigenfunction  as  a  side  benefit. 

Flammer  (1957)  and  Abramowitz  and  Steagun  (1965)  provide  power  series  approximations 
to  the  eigenvalue  that  are  valid  for  small  to  moderate  values  of  P,  the  dimensionless  electrical 
length.  These  approximations  provide  an  initial  input.  Flammer  related  the  eigenfunction  to 
an  infinite  series  of  associated  Legendre  polynomials.  A  three-term  recursion  relation  exists 
between  the  coefficients  of  the  series.  Starting  with  the  initial  approximation  and  the  term 
that  relates  the  angular  function  to  the  lowest  order  Legendre  polynomial,  we  recur  to  ever 
higher  order.  We  changed  the  eigenvalue  slightly  to  evaluate  the  derivative  of  the  highest 
order  coefficient.  We  then  changed  the  eigenvalue  to  minimize  the  value  of  the  last  term. 

(five  terms  would  have  been  sufficient  for  four-digit  accuracy.  The  1 5  terms  gave  us  every 
digit  listed  in  Flammer’s  book).  We  then  used  the  eigenvalue  and  coefficients  and  the 
normalization  condition  to  evaluate  the  angular  function. 

The  evaluation  of  the  derivative  of  the  angular  function  with  respect  to  r|  was  straight¬ 
forward.  We  used  the  expression  for  the  derivative  of  Legendre  polynomials  in  terms  of 
higher  order  ones.  We  checked  the  routine  to  an  accuracy  of  no  less  than  five  digits  by 
numerically  differentiating  the  expression  for  the  angular  function. 

The  expression  for  the  regular  part  of  the  radial  function  is  given  in  terms  of  spherical 
Bessel  functions  and  the  same  coefficients  found  for  the  angular  part.  Unfortunately,  as 
expressed  in  other  literature,  the  series  for  the  irregular  part  in  terms  of  spherical  Neumann 
functions  does  not  converge.  We  evaluated  the  integrals  given  on  pages  53  and  54  of 
Flammer’s  book  (equations  (50a)  and  (50b)  in  the  text)  by  Simpson’s  rule  at  one  value  of  u 
(usually  1.077).  The  results  for  the  real  and  imaginary  parts  of  the  function  and  its  derivative 
reproduced  every  digit  listed  in  Flammer’s  book.  Only  the  values  of  the  real  part  of  the  first 
four  modes  could  be  reproduced  to  six  digits  by  the  series  method.  We  do  not  feel  justified  to 
use  any  higher  modes. 

Using  the  values  of  the  integral  and  its  derivative  at  one  point  as  input,  we  solved  the 
differential  equation  for  the  radial  function  by  Runge  Kutta  methods.  We  compared  our 
analytic  solution  and  found  great  accuracy  for  almost  all  values  of  u. 

We  computed  the  integral  for  the  energy  used  in  calculating  the  value  of  Q  by  Simpson’s 
rule.  We  used  finite  values  of  u*  and  Uj+i.  We  then  compared  the  values  of  the  integral  as  the 
upper  limit,  Uj  became  large.  When  the  value  of  the  integral  did  not  change  in  the  digit,  we 
regarded  the  integral  as  having  converged  to  a  limit  for  an  infinite  domain. 
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